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Abstract 



We consider the Ujj^ —>■ oscillation measurements that would be possible at up- 
graded 1 GeV and multi-GeV conventional neutrino sources driven by future megawatt- 
scale proton drivers. If these neutrino superbeams are used together with detectors 
that are an order of magnitude larger than those presently foreseen, we find that the 
sensitivity to v^. oscillations can be improved by an order of magnitude beyond 

the next generation of accelerator based experiments. In addition, over a limited region 
of parameter space, the neutrino mass hierarchy can be determined with a multi-GeV 
long baseline beam. If the Large Mixing Angle MSW solution correctly describes the 
solar neutrino deficit, there is a small corner of allowed parameter space in which maxi- 
mal CP-violation in the lepton sector might be observable at a 1 GeV medium baseline 
experiment. Superbeams with massive detectors would therefore provide a useful tool 
en route to a neutrino factory, which would permit a further order of magnitude im- 
provement in sensitivity, together with a more comprehensive check of CP-violation 
and the oscillation framework. 
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1 Prologue 



We have recently completed a six-month study of the prospective physics program at 
a neutrino factory, evaluated as a function of the stored muon energy (up to 50 GeV) 
and the number of useful muon decays per year (in the range from 10^^ to 10^^). 
The basic conclusions presented in our report Q were that: (1) There is a compelling 
physics case for a neutrino factory with a muon beam energy of about 20 GeV or 
greater, and (2) The neutrino factory should provide at least 0(10^^) useful decays per 
year initially, and ultimately at least 0(10^*^) decays per year. The oscillation physics 
that can be pursued using initial electron-neutrino (z^e) and electron-antineutrino (Fg) 
beams provides the primary motivation for a neutrino factory. In particular we found 
that with 2 x lO^'^ decays per year, after a few years of running: 

(i) A — > oscillation signal could be observed, and the associated amplitude 
parameter sin^ 20i3 measured, for oscillation amplitudes approaching 10~^, three 
orders of magnitude below the currently excluded region and two orders of mag- 
nitude below the region expected to be probed by the next generation of long- 
baseline accelerator based experiments. 

(ii) Once a — > signal has been established in a long baseline experiment, matter 
effects can be exploited to determine the sign of the difference between the squares 
of the neutrino mass eigenstates 5777.32, and hence determine the neutrino mass 
hierarchy. 

(iii) If the large mixing angle MSW solution describes the solar neutrino deficit, and 
if sin'^ 2013 is not less than one to two orders of magnitude below the currently 
excluded region, a comparison of Ve i^^i and Ve oscillation probabilities 
would enable the measurement of (or stringent limits on) CP-violation in the 
lepton sector. 

(iv) Measurements of, or stringent limits on, all of the observable fe fx oscillation 
modes together with the observable vx modes, would enable a comprehen- 
sive test of the assumed oscillation framework. 

In parallel with our neutrino factory physics study, a companion design study 
was conducted to determine the feasibility of constructing a neutrino factory, and 
to identify the associated R&D issues. The design study concluded that a neutrino 
factory with the desired parameters was indeed feasible, although it would require a 
vigorous and well supported R&D program. Recognizing that the R&D would take 
some time, and that a neutrino factory would require a very intense (megawatt-scale) 
proton driver, it is reasonable to consider the neutrino oscillation physics program that 
could be conducted using a MW-scale proton driver en route to a neutrino factory. In 
our neutrino factory physics study report we recommended that an additional study 
of the oscillation physics potential at these "neutrino superbeams" be undertaken. 

The present document, which can be considered as an addendum to our initial 
report, presents results from a study of oscillation physics at 1 GeV and multi-GeV 
neutrino superbeams. 
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2 Introduction 



In this report we consider the oscillation physics capabilities of neutrino "superbeams" , 
which we define as conventional neutrino beams produced using megawatt-scale high- 
energy proton drivers. Examples of appropriate proton drivers are (i) the proposed 
0.77 MW 50 GeV proton synchrotron at the Japan Hadron Facility (JHF) Q, (ii) 
a 4 MW upgraded version of the JHF, (iii) a new ~ 1 MW 16 GeV proton driver 
that would replace the existing 8 GeV Booster at Fermilab, or (iv) a fourfold intensity 
upgrade of the 120 GeV Fermilab Main Injector (MI) beam (to 1.6 MW) that would 
become possible once the upgraded (16 GeV) Booster was operational. The 4 MW 
50 GeV JHF and the 16 GeV upgraded Fermilab Booster, are both suitable proton 
drivers for a neutrino factory. Hence a neutrino superbeam might provide a neutrino 
physics program en route to a neutrino factory. 

The next generation of accelerator based long-baseline neutrino oscillation experi- 
ments are expected to confirm the — > Vr oscillation interpretation of the atmospheric 
muon-neutrino deficit, and begin to measure the associated oscillation parameters with 
modest statistical precision. To make a significant improvement in oscillation measure- 
ments beyond the next generation of experiments will require a significant increase in 
signal statistics. A factor of a few increased neutrino beam fiux will not be sufficient 
unless there is also a substantial increase in detector mass. We will therefore assume 
that superbeam experiments will use detectors an order of magnitude larger than those 
currently under construction. Hence, a superbeam experiment would yield data sam- 
ples with a statistical sensitivity a factor of at least \/40 better than expected for 
K2K @, MINOS |], OPERA 0, and ICARUS [|. Clearly, this would enable signifi- 
cant progress in pinning down Uj- oscillations, for example. However, the neutrino 
oscillation physics program at a superbeam would have to justify the substantial in- 
vestment associated with the detector. For example, a detector costing ten times the 
MINOS detector would be of order $300M. Hence, more precise measurements of the 
quantities already measured by the next generation experiments may be insufficient 
motivation for a superbeam. 

The primary motivation for a neutrino superbeam is likely to be the search for, 
and measurement of, — > v^. oscillations. The observation of this mode would enable 
the associated amplitude parameter sin^ 2^13 to be determined, and open the way for 
the determination of the sign of the differences between the squares of the neutrino 
mass eigenstates ^^'^ hence the determination of the neutrino mass hierarchy. A 

comparison between the oscillation probabilities for — > and 17^ — > oscillations 
might also be sensitive to CP-violation in the lepton sector. Hence, in principle the 
i^fi Vf. mode at a superbeam can offer a handle on much of the physics that the 
T^ii mode offers at a neutrino factory. Indeed, for low energy neutrino beams 
{Ey < 10 GeV) the neutrino fiuxes at a superbeam are comparable or larger than the 
corresponding neutrino factory fluxes (see Table 1). However, there is a substantial 
qualitative difference between searching for i/g appearance at a superbeam and i^^ ap- 
pearance at a neutrino factory. The signal signature at a superbeam is the appearance 
of an isolated electron (or positron) in a charged current (CC) event. As we will see, 
this signature is plagued with backgrounds at the level of 0(1%) of the total CC rate. 
For comparison, the signal signature at a neutrino factory is the appearance in a CC 
event of a wrong-sign muon (a muon of the opposite charge-sign to that of the muons 
stored in the muon ring). This signature enables backgrounds to be suppressed to the 



4 



Table 1: Neutrino event rates assuming no oscillations, compared with intrinsic beam back- 
grounds for conventional and muon-derived beams of comparable energies. The calculations 
assume a 1.6 MW proton source is used for the MINOS-type beam, the neutrino factories 
provide 2 x 10^° muon decays per year in the beam-forming straight section, and the detector 
is 732 km downstream of the neutrino source. 



Beam 


<E^> 


CC Events 




(Signal: i/^ Ve) 


(GeV) 


(per kton-year) 


Fraction 


MINOS-LE 


3.5 


1800 


0.012 


MINOS-ME 


7 


5760 


0.009 


MINOS-HE 


15 


12800 


0.006 




Beam 


<E^> 


Ve CC Events 




(Signal: v^) 


(GeV) 


(per kton-year) 


Fraction 


4.5 GeV /U Ring 


3.5 


400 





9.1 GeV ^ Ring 


7 


3700 





18.2 GeV Ai Ring 


15 


31400 





30 GeV iJL Ring 


20 


72600 






level of 0(0.01%) of the total CC rate Q. Hence, to understand the oscillation physics 
potential at a superbeam we must have a good understanding of the backgrounds and 
the systematic uncertainties associated with the background subtraction. 

In Section 3 of this report we begin by discussing the properties of conventional 
neutrino beams. The motivation for large mass detectors with excellent background 
rejection is discussed in Section 4. The most important backgrounds to i/^ — > Ue and 
z?^ — > i>f> oscillations are discussed in Section 5. The physics capabilities of multi- 
GeV long baseline experiments and 1 GeV medium baseline experiments are discussed 
respectively in Sections 6 and 7. A summary is given in Section 8. Throughout we 
will use the three-flavor oscillation framework which is reviewed in Appendix 1, with 
the "leading" oscillation parameters sin^ 2^23 and 5m\2 determined by the atmospheric 
neutrino deficit, and the sub-leading parameters sin^ 29 12 and dm^i determined by the 
solar neutrino deficit. This choice is appropriate if the LSND effect Q is not confirmed 
(by MiniBooNE pO|] , for example). Should the LSND effect be confirmed, there is 
likely to be a strong physics case for a low intensity neutrino factory, which might be 
constructed on a relatively fast timescale with only a short R&D phase. The case for 
a separate superbeam program is less obvious (or at least different) in this case. 

3 Beam characteristics and event rates 

A conventional neutrino beam is produced using a primary proton beam to create a 
secondary beam of charged pions and kaons, which are then allowed to decay to pro- 
duce a tertiary neutrino beam. The secondary beam can be charge-sign selected to 
produce either a neutrino beam from positive meson decays or an antineutrino beam 
from negative meson decays. The secondary particles, which are confined radially us- 
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Figure 1: Differential spectra for pions (left) and kaons (right) produced when 120 GeV 
protons are incident on the MINOS target. The distributions are normalized to correspond 
to the number of particles per incident proton. 

ing either a quadrupole channel or horn focusing, are allowed to decay in a long decay 
channel. The resulting neutrino beam consists mostly of muon neutrinos (or antineu- 
trinos) from vr^ — > //f^ decays, with a small "contamination" of electron neutrinos, 
electron antineutrinos, and muon antineutrinos from muon, kaon, and charmed meson 
decays. The fractions of Vg, T^e and in the beam depend critically on the beamline 
design. 

Figure |l] shows the momentum spectrum of charged pions and kaons produced when 
a 120 GeV beam of protons strikes a 2 interaction length graphite target [11|. Note that 
at low momenta, positive and negative secondaries are created at comparable rates, 
but at higher momentum there is a marked asymmetry. For example, at 20 GeV the 
ratio of positive to negative secondaries is 3/2. Hence, for high energy beams there is 
a flux penalty in producing an antineutrino beam rather than a neutrino beam. This 
flux penalty increases with increasing beam energy. 

For the two-body decays vr — > fii'^ and K fiv^ there is a one to one correspon- 
dence between the energy of the parent meson and the energy of the neutrino at the 
far detector. For a parent particle of mass nih and energy Eh, traveling in a direction 
9i,fi with respect to the far detector, the neutrino energy is given by: 



mr, 



rrih 



2mh 



Eh - Ph cos e 



Eh 



mr 



m. 



uh 



mt 



1 + 7 



2fl2 

vh 



(1) 



For a perfectly focused beam 9yh = 0, and = 0A2Eh for pion decays and 0.95Eh 
for kaon decays. In practice the beamline is designed to focus pions within a given 
momentum window. A broader pion momentum acceptance will result in a higher flux 
of neutrinos in the forward direction, but will also result in a broader energy spread 
within the neutrino beam. The flux of neutrinos per meson at the far detector is given 
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by: 



where BR is the branching fraction for the appropriate meson decay, 7 is the Lorentz 
boost of the decaying particle, and L is the distance to the detector. Note that the 
flux at the far detector has the famihar 7^ dependence, and since the cross-section 
increases hnearly with 7, the event rate has a 7^ dependence. Figure ^ shows, for a 
perfectly focussed beam, the calculated and i/^ event rates per kton-year for the 
MINOS detector at L = 732 km, assuming a total decay region of 725 m. Note that 
for a realistic focusing system the neutrino flux would be reduced, typically by a factor 
of 2 or 3. The calculation shown in Fig. ^ is for 15 x 10^'' 120 GeV protons striking 
a graphite target (4 x NuMI [12| for 1 year) similar to the one being constructed for 
the MINOS beam |11], with the beamline set to focus either positive or negative pions. 
Also shown are the corresponding event rates at a neutrino factory in which there are 
2 X 10^0 useful muon decays, with i?^ = 10 and 20 GeV. 

The total CC event rates for 1.6 MW NuMI low, medium, and high energy su- 
perbeams are compared in Table |l| with the corresponding (same < Ep >) rates at a 
neutrino factory. Note that the rates at a neutrino factory rapidly exceed the corre- 
sponding conventional beam rates for neutrino beam energies exceeding about 10 GeV. 
For lower energies, conventional beams provide higher rates although the beam is not 
background free, and the electron appearance signal is experimentally challenging. 

Finally, the decay channel for a very long-baseline superbeam must fit within the 
viable rock layer below the proton driver. At Fermilab this rock layer is ~ 200 m deep, 
below which there is a deep aquifer. This limits the length of the decay pipe to less 
than 200m X sin 0, where the dip angle 9 depends on the baseline: L = 12756 x sin 9 km. 
In fact the length of the decay region is further restricted by the depth of rock used to 
bury the proton driver, bend the proton beam to the required direction, accommodate 
the target and focusing systems, and if there is a near detector, accommodate the 
associated shielding and near detector hall. Figure ^ shows the relative flux loss for 
different energy beamlines as a function of baseline length |13|. The calculation allows 
30 m for the near detector shielding plus hall. For a baseline of 732 km, 9 = 3.3° 
and the NuMI beam pipe length of 675 m is not restricted by the depth of the good 
rock. However for a far site at 7300 km (Fermilab Gran Sasso), 9 = 35°, and 
the length of the decay channel is severely restricted, so that for the LE (ME) [HE] 
NuMI beams only 50% (25%) [13%] of the pions would decay in the channel. Clearly, 
decay channel length restrictions must be taken into account in calculating the fluxes 
for very long baselines (> 3000 km). For trans- Atlantic or trans-Pacific baselines 
there is a premium on minimizing the shielding and detector hall length for the near 
detector, minimizing the length of the targeting hall, using high-field dipoles to rapidly 
bend the 120 GeV proton beam to the required direction, and perhaps considering a 
"roller-coaster" geometry for the proton beam. 



4 Large detectors and low backgrounds 

Our first physics goal is the observation of — > Ve oscillations. Since, to a good 
approximation the oscillation probability is proportional to the amplitude parameter 
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Figure 2: Differential z/^(top) and (bottom) event rates in the MINOS detector 732 km 
downstream of a perfectly focused beam of pions and kaons produced when 15 x 10^° 120 GeV 
protons are incident on a graphite target (4x NuMI for 1 year). The top panel also shows 
the predicted spectrum for a realistic focussing system, namely the MINOS HE beam (as 
indicated). Note that the rates arc expected to be a factor of 2 or 3 less for a realistic system 
than for a perfectly focussed beam. The distributions are compared to the corresponding 
i^e event rates 732 km downstream of 10 GeV and 20 GeV muon storage rings driven by a 
1.6 MW 16 GeV proton source. 
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2000 4000 6000 8000 10000 
Baseline (km) 

Figure 3: Fraction of tt^ that decay in a channel of maximum length within a 200 m deep 
rock layer, shown as a function of baseline for three different average NuMI beam energies: 
3.5 GeV (LE), 7 GeV (ME), and 15 GeV (HE). For convenience, distances from Fermilab to 
SLAC/LBNL, Gran Sasso, and Kamioka are indicated by vertical broken lines, and the decay 
fractions for the foreseen NuMI beams {L = 730 km) are indicated by dotted horizontal lines. 

sin^ 2^13, it is useful to define the sin^ 20i3-reach for a given experiment, which we define 
as that value of sin^ 2^13 which would result in a signal that is 3 standard deviations 
above the background. Taking the atmospheric neutrino deficit oscillation scale 6711^2 
to be in the center of the region indicated by the SuperKamiokande (SuperK) data, for 
a given proton driver, superbeam design, and baseline, we can calculate the sin^ 2^13 
reach once we specify (a) the data sample size D (kt-years), defined as the product 
of the detector fiducial mass, the efficiency of the signal selection requirements, and 
the number of years of data taking, (b) the background fraction fs, defined as the 
background rate divided by the total CC rate for events that pass the signal selection 
requirements, and (c) the fractional uncertainty c/^/Zs. Note that D determines 
the statistical uncertainty on the signal. For a fixed D, fs determines the statistical 
uncertainty on the background, and Jb x <^fB/ Ib determines the systematic uncertainty 
on the background subtraction. 

Contours of constant sin^ 2^13 reach that correspond to various values of 
are shown in the (/b, L>)-plane in Fig. ^ for Ey = 1 GeV superbeams produced using 
0.77 MW and 4 MW JHF proton drivers. The corresponding contours are shown in 
Fig. ^ for long-baseline 1.6 MW medium energy and high energy NuMI superbeams. 
The contours have a characteristic shape. At sufficiently large D the sin^ 2^13 sensitivity 
is limited by the systematic uncertainty on the background subtraction, and the reach 
does not significantly improve with increasing dataset size. The contours are therefore 
vertical in this region of the figures. At sufficiently small D the sensitivity of the 
i^fj, appearance search is limited by signal statistics, and further reductions in 

Jb do not improve the sin^ 2^13 reach. The contours are therefore horizontal in this 
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Figure 4: Contours of constant sin^ 26*13 reach that correspond to a z/g — > t'^ signal that is 
3 standard deviations above the background The contours are shown in the (D,/^)- 
plane, where D is the data-sample size and the background rate divided by the total CC 
rate. The contours are shown for the 0.77 MW (left-hand plots) and 4.0 MW (right-hand 
plots) JHF scenarios with L = 295 km. The top panels show curves for (Jf^/ fs = 0.1, while 
the bottom panels show curves for cffg/ fs = 0.1, 0.05, and 0.02. 
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Figure 5: Contours of constant sin^26'i3 reach |13] that correspond to a z/g ^ z/^ signal that 
is 3 standard deviations above the background, at upgraded 1.6 MW NuMI ME (left) and 
HE (right) beams. The contours are shown in the (D, /B)-plane, ■where D is the data-sample 
size and the background rate divided by the total CC rate. The contours are sho"wn for 
L = 2900 (top), 4000 (center), and 7300 km (bottom). Curves are sho^wn for systematic 
uncertainties on the background subtraction <Jfg/ fs = 0.1, 0.05, and 0.02. 
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region of the {/b, D)-plane. Note that the next generation of neutrino experiments are 
expected to achieve sin^ 2^13 reaches of ~ 0.03. Now consider as examples the 4 MW 
JHF beam, and the 1.6 MW NuMI beam with L = 2900 km. Taking ct/^/Zb = 0.05, 
an inspection of the figures leads us to conclude that if we wish to improve the sin^ 2^13 
sensitivity an order of magnitude beyond the next generation of experiments, in the 
limit of very massive detectors at the JHF (NuMI) superbeams, we can only tolerate 
background fractions Jb < 0.007 (0.004). Achieving the required background rejection 
in a water cerenkov detector seems challenging. Achieving the required detector masses 
with detector technologies that can meet the /b requirements also seems challenging. 
Hence an understanding of the physics capabilities at a superbeam necessarily begins 
with an understanding of the parameters D, /b, and CTfg/fB that can be achieved with 
realistic (but futuristic) detectors. 

In the following section we consider the backgrounds, and for realistic futuristic 
detectors, D and Jb- This will enable us to use the {fB,D) figures to determine the 
sin^ 2^13 sensitivity at various superbeams. 

5 Backgrounds to z/^ z/g oscillations 

Backgrounds play a critical role in determining the sensitivity of a superbeam experi- 
ment to — > oscillations. In our neutrino factory studies ||l| it was straightforward 
to obtain backgrounds to the Ue — > i^^ search at the Jb ~ 10~^ level by cutting on 
the final state muon momentum. However, for conventional neutrino beams reduc- 
ing the background rates below 10~^ of the total CC rate is not trivial and would 
require significant reductions in the beam flux and signal selection efficiency. There 
are four important sources of background in a i/g appearance search at a superbeam: 
(i) electron-neutrinos produced in the initial beam, (ii) neutral current (NC) neutrino 
interactions in which a tt^ is mis-identified as an electron, (iii) CC interactions in 
which a tt^ is mis-identified as an electron, and the muon is not identified, and (iv) 
events from Vr oscillations followed by v^- CC interactions and either t ^ e + X 

decays, or r — > tt^+X decays in which the vr'' fakes an electron. In the following we dis- 
cuss these backgrounds and how, for different detector technologies, the backgrounds 
might be suppressed. 

5.1 Electron— neutrino contamination 

In a conventional neutrino beamline, muon neutrinos are produced in two-body pion- 
and kaon-decays. However, the charged and neutral kaons can also undergo 3-body 
decays to produce an electron neutrino (or antineutrino) . Furthermore, secondary 
muons in the beamline can also decay to produce electron neutrinos. Generally speak- 
ing, higher primary proton beam energies yield higher kaon/pion ratios, and longer 
decay channels permit more muon decays. The expected electron neutrino background 
fractions are listed in Table || for the next generation of conventional neutrino beams. 
Note that for most of the beamlines listed the background fraction is around the 
1% level. The exceptions are the MiniBooNE beam (which benefits from a relatively 
low primary proton energ y), and ORLaND [|l| (which uses stopped muons). Since the 
electron neutrinos within the beam are created in three-body decays, the energy 
spectrum is typically much broader than the spectrum. This effectively means that 
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Table 2: Electron neutrino fractions and the fractional energy spreads for a selection of 
current (or next) generation conventional neutrino beams. Note that most beamlines produce 
a beam with a fractional energy spread between 30% and 50%, and a v^. contamination that 
ranges from 0.2% to 1.2%, for beams at or above 1 GeV. 



Beamline 


Proton 
Energy (GeV) 


Peak 
Energy (GeV) 


ratio 




K2K 


12 


1.4 


0.7% 


1.0 


MINOS LE 


120 


3.5 


1.2 


% 


0.28 


MINOS ME 


120 


7 


0.9 


% 


0.43 


MINOS HE 


120 


15 


0.6 


% 


0.47 


CNGS 


400 


18 


0.8 


% 


0.33 


JHF wide 


50 


1 


0.7 


% 


1.0 


JHF HE 


50 


5 


0.9 


% 


0.40 


MiniBooNE 


8 


0.5 


0.2 


% 


0.50 


ORLaND 


1.3 


0.0528 


0.05% 


0.38 



there is always some fraction of the electron neutrino flux which overlaps in energy 
with the muon neutrino flux. 

The fe flux contributions from and /i^ decays can be reduced by decreasing the 
secondary particle momentum acceptance. The contributions from decays can be 
reduced by putting large bends in the beamline immediately after the proton target. 
The neutral meson backgrounds are particularly dangerous for v running because they 
give electron neutrinos in a Ve appearance search, and the cross section for the 
background is therefore twice as large as for the Ue signal. It has been suggested [15| 
that a neutrino beam could be made with an extremely small contamination by 
only accepting pions and kaons within a narrow momentum interval, and then rejecting 
neutrino events that have a total energy inconsistent with the expected neutrino beam 
energy. To understand how this might work we consider the correspondence between 
the the neutrino beam energy spread and the contribution to Jb from the initial 
flux, which comes from K^, /j,^, and decays. The ratio of the initial flux to 
flux is shown as a function of the fractional beam energy spread in Fig. |6| for a perfectly 
focussed secondary beam produced with 120 GeV primary protons on a 2 interaction 
length graphite target, followed by a 725 m long decay channel |11]. To achieve a 
background fraction that is no larger than 0.1% in either u or i> running would require 
a beamline with a momentum acceptance no larger than 10%. In addition, we must 
suppress the remaining z^g contribution from decays (using one or more dipoles after 
the proton target) by factors in excess of 2.5 and 5 for v and i? running respectively. 
Noting that the effective flux also decreases roughly linearly with the momentum 
acceptance, we conclude that, even for an idealized perfectly focussed beam, reducing 
the initial z/g contamination in the beam to 0.1% will result in an order of magnitude 
reduction in D. 

An alternative way of decreasing the contamination from kaon decays is to de- 
crease the primary proton energy (from 120 GeV to something less). However, the 
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Figure 6: Fraction of z^e + Pe events in a (left) or z/^ (right) beam shown as a function of 
the fractional beam momentum spread. The various symbols correspond to different mean 
beam energies, from 5 GeV (filled circles) to 17 GeV (open squares). The distributions that 
increase (decrease) from left to right are the contributions from K^,fi^ {Kl) decays. 

MiniBooNE beam, which uses 8 GeV primary protons, is expected to achieve a 
contamination of no better than 0.2%. Furthermore, even if kaon contributions are 
completely ehminated, muon decays still provide electron neutrinos. 

In summary, for a multi-GeV neutrino beam, we would not expect a reduction in the 
momentum acceptance of the decay channel to significantly reduce the contamination 
in the beam unless we are willing to accept a large reduction in D. Optimistically, 
the contribution to Jb from the Ve contamination in the beam might be reduced to 
~ 0.2 — 0.5% for the multi-GeV neutrino beams considered later in this report. 



5.2 Neutral Current Backgrounds 

For high energy neutrinos the NC cross section is roughly 40% of the CC cross section, 
and is independent of the initial flavor of the interacting neutrino. A NC interaction 
can result in a total visible energy anywhere from zero to the initial neutrino energy 
Ey. Some of the visible energy may be in the form of neutral pions, which can fake 
a single electron signature. The probability for a NC event to produce an energetic 
vr'^ that fakes an electron depends both on the vr'^ production rate and on the details 
of the detector and signal selection requirements. In the following we will begin by 
considering tt^ production in NC events, and then consider various detector strategies 
for reducing the background. 



5.2.1 Energetic 7r° production 

Consider the probability that a NC interaction will produce an event with an energetic 
vr'^ that could fake a Ve CC interaction. Our calculations have been performed using 
the GEANT Monte Carlo program together with a LEPTO simulation to generate 
neutrino interactions ll^. A rough approximation has been used for quasi-elastic and 
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Figure 7: Energy distributions for electrons and neutral pions produced in 25 GeV neutrino 
interactions. Solid line: Electrons produced in CC interactions (scaled by 0.04). Dotted line: 
Neutral pions generated in NC interactions. Dashed line: Electrons from charm semileptonic 
decay. 

resonance production in CC events, but no resonance production has been included 
for the NC events. Note that the lower the neutrino energy, the higher the fraction of 
quasi-elastic and resonance interactions. As a benchmark, for a 5 GeV neutrino beam 
the quasi-elastic and resonance contributions are about 33% of the total event rate. 

The simulated energy distribution for electrons produced in 25 GeV Ue CC events 
is compared in Fig. ^ with the corresponding distribution for vr'^'s produced in 25 GeV 
NC events. The -k^ rates are small at high energy. A cut on the energy of the electron 
candidate can therefore reduce the NC background in a ^ Ve search. The most 
dangerous background events are those in which the vr'^ takes most of the energy of the 
hadronic fragments, and hence x = E^o /£'had is large. These neutral pions can fake 
an isolated electron from a CC interaction. The fraction of NC events having a vr'^ 
with energy greater than a given fraction of i?had is shown in Fig. ^ for different ranges 
of hadronic energy. The vr'^ fragmentation functions are roughly independent of -Ehad- 
To a good approximation, at large x 0.3) a single function describes all the curves 
shown in the figure: 

p{x) = (0.49) - (0.96)2; + (0.47)^2 (3) 

Note that in a NC event -Ehad is just the visible energy. Within the framework of the 
parton model, for a given neutrino energy, the visible energy spectrum is described 
by HI: 

1 dN 15 ( ^ {^-yf\ 



N dy 16 V ^ 



(4) 



where y = £^had/^v 

We can use Eqs. |^ and Q to calculate the probability that there will be a vr'' above 
some cut-off energy. Thus, for a given input neutrino energy spectrum, an output NC 
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Figure 8: Fraction of NC events with a rc^ with energy greater than a given fraction of 
the hadronic energy. The different histograms are for different ranges of hadronic energy: 
Sohd: 20 < E^^ < 25 GeV. Dashed: 15 < Ehad < 20 GeV. Dotted: 10 < Ehad < 15 GeV. 
Dot-Dashed: 5 < -Ehad < 10 GeV. The curve shows the parameterization given in the text. 

background spectrum can be calculated under the assumption that all of the vr'^'s are 
incorrectly identified as electrons. As an explicit example, consider a two-horn (NuMI- 
like) neutrino beam with a mean energy of 10 GeV, and a baseline of 9300 km (Fermilab 
— > SuperK), and assume the oscillation amplitude parameter sin^ 2^13 = 0.01, and the 
leading oscillation scale 5m'^2 — 3-5 x 10~^ eV^. The expected V(, signal is 

compared in Fig. ^ (left panel) to various background components. In particular, the 
NC background spectrum is shown after requiring that the electron-candidate energy 
is at least 30% of the total energy, which suppresses the vr'^'s. Note that the surviving 
NC background is still the dominant background component, although the others are 
not negligible. The total background level is about equal to the signal level for this 
particular example. The signal/background ratio could be improved by narrowing the 
energy spread of the neutrino beam. For example, if neutrinos with energies between 
8 GeV and 10 GeV are selected, the right-hand panel in Fig. ^ shows the resulting 
signal and background distributions. Note that for this case the NC background level 
is roughly 1/6 of the signal. 

In summary, by requiring the electron candidate to have a significant fraction of 
the visible event energy, the NC background can be reduced so that its contribution 
to fs ~ 0.04, at the price of reducing D by ~ 30%. Further reductions in the NC 
background depend on detector details, and will be discussed in the following sub- 
sections. 

5.2.2 Liquid Argon TPC 

In a liquid argon TPC there are four ways that NC events can contribute backgrounds 
in an electron appearance search: a) a photon converts very close to the interaction 
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Figure 9: Differential distributions for NC, Ut-, and intrinsic z/g backgrounds after requiring 
the electron candidate has at least 30% of the total event energy, compared with corre- 



sponding distributions for a 



signal (sin^2^i3 = 0.01, 



3.5 X 10-3 eV^ 



at L = 9300 km. In the left panel a NuMI-like beam with a mean energy of 10 GeV is 
assumed. The right panel shows only the contributions from neutrinos in the beam with 
energies between 8 GeV and 10 GeV. Solid line: CC signal. Dashed line: NC background. 
Dot-Dashed line: Ur background. Dotted line: Intrinsic Ue background. 

vertex, b) a charged pion interacts very close to the interaction vertex, c) a charged pion 
overlaps with a vr'' creating an electron-like track, or d) there is an asymmetric Dalitz 
decay or external conversion of a vr'' created at the interaction vertex. The imaging 
capability of a hquid argon TPC of the type planned by the ICARUS coUaboration 
is expected to facihtate much better rejection of these backgrounds than any other 
high-mass neutrino detector that exists or is planned. 

A full simulation of the ICARUS detector g in the CNGS beam (which has a 
mean neutrino energy of about 17 GeV) shows that background sources a) and d) have 
the largest rates, and contribute respectively 4262 and 275 out of 15550 NC events. 
Note that there are two charged particles resulting from the Dalitz decay or the photon 
These backgrounds can therefore be suppressed by requiring that the 



conversion. 



electron candidate does not have an unusually high dE/dx before it showers. Figure 10 



shows an ICARUS Monte Carlo simulation of the electron/Dalitz pair separation using 
the first 6 cm of the "track" after the interaction vertex. Requiring the deposited 
energy to be less than 1 MeV would keep 90% of the signal while removing over 99% 
of the two-electron background. Requiring the reconstructed electron energy to exceed 
1 GeV would further reduce these backgrounds by about a factor of 3. Hence, before 
any other kinematic or reconstruction cuts are applied, backgrounds a) and d) can be 
reduced to ~ 10^^ of the NC event rate 

The NC background sources b) and c), where a single charged pion either interacts 
or overlaps with a vr^, have a low rate (31.5/15550 NC events). Once again, two thirds 
of these events can be removed by requiring that the electron candidate has an energy 
greater than 1 GeV. Hence, backgrounds b) and c) can be reduced to less than 10"'^ 
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ure 10: Simulated dE/dx measurement in liquid argon for both single electrons and Dalitz 
based on summing over 20 wires in the ICARUS detector, or 6 cm of Liquid Argon. 

of the NC event rate, at the cost of reducing D by a factor of ~ 0.9. 

We conclude that, in an ICARUS-type detector, the total contribution to the back- 
ground fraction fs from NC events can be reduced to ~ 0.001, at the cost of a modest 
reduction in the dataset size Z) of ~ 10 — 15%. 



5.2.3 Water Cerenkov 

Although Super-Kamiokande can easily discriminate between quasi-elastic Ve and quasi- 
elastic events, this becomes problematic in the presence of a hadronic shower with 
a significant number of high energy secondaries. It was shown Q that for neutrino 
energies of ~ 1 GeV, the NC background contribution to fs can be reduced to ~ 0.03 
at the cost of a reduction in D by a factor of 0.68. We will assume that, with further 
optimization, Jb can be reduced to ~ 0.02. The vr'^'s in NC events are rejected by 
exploiting the cerenkov ring structure (two separated rings from vr*^ 77 decays, one 
ring for an electron shower). However, above neutrino energies of ~ 1 GeV the two 
photons from the decay are indistinguishable, and the only discrimination against 
NC backgrounds is based on the electromagnetic energy fraction in the event |p!s|]. 
Hence we conclude that for 1 GeV neutrino beams the contribution to /b from NC 
events can be reduced to about 0.02. For higher energy beams the backgrounds are 
much larger, and we will assume that /b = 0.04. 



5.2.4 Sampling Calorimeters 

Sampling calorimeter neutrino detectors typically consist of a sandwich of target ma- 
terial (Iron, glass. Lead) and active material (scintillator or wire chambers). Relatively 
crude sampling calorimeters offer a compact and cost effective and way of instrument- 
ing a large fiducial mass for muon detection, and hence CC measurements. However, 
a — > z^e experiment requires an electron appearance measurement, which imposes 
significant additional requirements on the granularity of the calorimeter. In order 
to distinguish the electron produced in a neutrino interaction from production. 
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the calorimeter must be finely segmented both longitudinally and transversely. The 
longitudinal segmentation is needed to resolve the difference between 1 particle (an 
electron) and 2 particles (from a photon conversion) at the beginning of the electro- 
magnetic shower. The transverse segmentation is needed to prevent accidental overlaps 
of charged particles which would veto a real electron. We consider both fine-grained 
and course-grained sampling calorimeters: 

(i) Fine grained calorimeter. The Soudan 2 detector is an example of a finely seg- 

mented iron calorimeter with drift tube readout. The steel plates are 3.2 mm 
(0.18 Xq) thick. The estimated fraction of NC events which fake fg interactions 
is 2.3% after applying cuts that include requiring the candidate electron energy 
to be at least 50% of the visible energy. Hence, the contribution to /b from NC 
events in a Soudan 2 type detector is ~ 0.01. The associated reduction in D is 
by a factor of 0.6. 

(ii) Course grained calorimeter. The MINOS detector is an example of a course 
grained iron calorimeter with scintillator readout. The toroidally magnetized 
steel plates are 1 inch thick, and the detector is instrumented with 4.1 cm wide 
scintillator strips. NC background events can be rejected based on the electron 
candidate energy, the total event energy, and the number of scintillator strips 
associated with the electron shower. A recent study |2C] suggests that with the 
medium energy NuMI beam, the NC backgrounds for a ^ i^e search in a 
MINOS-like detector can be reduced so that the contribution to fs is ~ 0.01, at 
the cost of reducing D by a factor of 0.33. 



5.3 Charged Current Backgrounds 

In general, in a f ^ — > Ue search most CC events can be rejected due to the presence 
of a muon. In our present study we only consider the vr^ backgrounds from NC events. 
However, for the NuMI ME beam our simulations show that if all muons below 2 GeV 
in CC events are missed, the vr*^ background from CC events will roughly be equal to the 
corresponding background from NC events. The muon detection threshold is detector 
dependent, and more detailed detector studies are needed for the CC background to 
be taken into account. 



5.4 Backgrounds from Ur oscillations 

The oscillation process Uj- can produce a i/^ — > i/g background if the Uj- interacts to 
produce a tau-lepton that subsequently decays electronically {t ^ e + X, BR = 20%). 
If the detector does not have good discrimination between electrons and neutral pions, 
the process r — > utt^Xi't- {BR = 37%) also contributes to the background. The Ur 
backgrounds are particularly dangerous because they have the same dependence on 
5m\2L/E as the signal. Figure |ll| shows the CC cross section ratio as a function 

of neutrino energy |21]. Clearly the v-j- backgrounds can be eliminated by running 
below, or near to, the CC threshold (5 GeV). 

If the neutrino beam energy is significantly above 5 GeV, one way to reduce the 
Vr background is to require that the electron candidate carries a significant fraction of 
the total observed interaction energy. Assuming a perfect detector energy resolution. 



the distributions of electron energy divided by total energy are shown in Fig. 12 for 
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Figure 12: Ratio of electron energy to total visible energy for (left) v^. and v^. charged current 
events and (right) and Vr charged current events, where the r subsequently decays to an 
electron. 
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Figure 13: For oscillations with (5m|2 = 3.5 x 10"'^ eV^ and sin^ 2^23 = 1, the i'^ — >■ z/g 
background contribution /b(t) arising u^- interactions followed by r — > e decays (Left panel), 
and the associated z/^ signal efficiency (right panel) shown as a function of a cut on 

the minimum electron energy. The energy resolution assumed for electrons is a{Ee)/Ee = 



5%/^E,{GeV). 

10 GeV i^e and CC events (left), and 10 GeV Vr and v'r CC events (right) where 
the r subsequently decays to an electron. If the electron candidates are required 
to carry at least 50% of the visible energy, the remaining i/j- background level (for 
sin'^(1.27(5m32-Z^/-E') = 1) is a few per cent at 10 GeV. Note that the backgrounds 
for and v' running are different. Better background rejection (for a given reduction 
in signal efficiency) can be obtained by imposing a minimum energy requirement on 
the electron candidate instead of a requirement on the fraction of the total energy 
carried by the electron. The contribution to /b from the r — > e + X backgrounds 
is shown in Fig. |l^ as a function of the minimum electron energy requirement (left 
panel) for the NuMI high energy beam with L = 7300 km, and for the medium energy 
beam with L = 2800 km. The corresponding reductions in D are shown in the right 



panel of Fig. 13. Table ^ summarizes the expected contribution to /b for three NuMI 
beams using some explicit electron energy cuts. For the medium- and high-energy 
beams, to reduce the contribution to /b to ^ 0.01 requires a significant reduction in 
D. Furthermore, the total contribution to /b from ly-r interactions can be a factor of ~ 2 
larger than listed in Table ^ if the detector does not provide adequate discrimination 
against r — > mr^Xur decays. 

There are other kinematic handles that can in principle be used to suppress i^r 
backgrounds, but the performance of a set of given kinematic cuts depends upon the 
detector details. Obvious kinematic quantities that can be exploited are the missing 
transverse momentum (penalizing the missing energy associated with the neutrinos 
produced in the r-lepton decay), and the electron transverse momentum distribution. 



The simulated distributions for these variables are shown |22] in Fig. 14 for several event 



types in the ICARUS detector at the CNGS beam. Note that the CNGS beamline has 
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Table 3: Summary of v^- backgrounds in which r — > e + X. The background rates are hsted 
for NuMI-Uke low-, medium-, and high-energy beams, assuming the oscillation parameters 
|(5m|2| = 3.5 X 10~^ey^, and sin^ 2^23 = 1- The assumed electron energy resolution is 
a{Ee)/Ee = 5%/ ^Ee{GeV), as anticipated for the ICARUS detector. Note that the total 
i>T- related backgrounds can be a factor of ~ 2 larger if the detector does not have good 
descrimination against r nvr^Xz/^ decays. 



Beam 


Distance 


fB{r) 


fsir) 


Ee 


Signal 




(km) 


no cut 


with cut 


cut 


Efficiency 


NuMI-LE 


732 


0.001 


0.001 


none 


1 


NuMI-ME 


2800 


0.013 


0.006 


> 3 GeV 


0.6 


NuMI-HE 


7300 


0.031 


0.013 


> 5 GeV 


0.6 



Ain = 3.5 X 10'^ eV^ sin 2&2^=0.9 sm20,,=O.l hni = 3.5 x 10'^ eV' sin 2® ^^=0.9 sin2@j,=0.1 
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Figure 14: Missing transverse momentum (left) and electron transverse momentum (right) 
distributions for z/e charged current events and for z/^, r — >• e events in the ICARUS detector, 
assuming oscillation parameters as indicated. 
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a Ve background fraction of 0.8% and a mean neutrino energy of 17 GeV, and sin^ 2^13 
has been assumed to be 0.1, just at the current hmit. The distributions shown in the 
figure are for all electron-like events with no electron energy requirement. It should be 
noted that events that pass an electron energy cut will tend to have electrons moving 
close to the Vr direction, and hence are more likely to survive other kinematic cuts. For 
an experiment in which the r appearance rate is likely to be significant, the ICARUS 
collaboration finds the best sensitivity is obtained with a combined fit, rather than by 
eliminating the Vr events. 

We conclude that backgrounds are only significant for the medium- and high- 
energy beams, in which case an electron energy cut and/or good r ^ e rejection is 
needed to reduce the contribution to fs to ^ 0.01. This might be accomplished in a 
liquid argon detector, which provides good discrimination against r mr^Xur decays. 
The r related backgrounds might be larger in other detectors. For example, a recent 



study Q concluded that, with = 0.003 eV^ and L = 2900 km, after cuts which 
reduced D by a factor of 0.33, the contribution to /b from Uj- associated inter- 

actions was 0.02 for a MINOS-like iron-scintillator detector at a NuMI-like medium 
energy beam. 

5.5 The systematic uncertainty on the background 

Since a — > t'e search will not be background free, the background must be subtracted 
from any potential signal. The background subtraction introduces a systematic un- 
certainty associated with the imperfect knowledge of the expected background rate. 
Assigning a systematic error on the background rate requires understanding the uncer- 
tainty on the Ve contamination in the beamline, and the uncertainty on the detector's 
ability to reject background events. 

To understand the fg fraction in the initial neutrino beam, we must know the 
charged and neutral kaon components in the secondary beam. This requires knowledge 
of their production at the proton target, and knowledge of the beamline acceptance. 
In addition, the backgrounds from muons decaying in the decay tunnel must be under- 
stood, although generally speaking the uncertainty from this background component 
is relatively small. Measuring the fraction of charged and neutral kaons produced in 
the target is the subject of much current study. The experiment P907 [^] is being 
proposed at Fermilab to measure the kaon production cross sections for 120 GeV pro- 



tons on the MINOS target. A similar experiment, HARP |24], is planned at CERN to 
study meson production for the CNGS beam. P907 and HARP are expected to sig- 
nificantly reduce the kaon production uncertainties for neutrino beams using 120 GeV 
and multi-GeV (up to 15 GeV) proton beams, together with targets similar to NuMI 
and CNGS. However, the combined effects of both production and acceptance must 
be understood. For this purpose, a very fine-grained near detector could be used. As 



an example, the NOMAD detector |25] can determine the neutral kaon contribution to 
the Ve flux by comparing the rates of Ve and events that are seen. With an electron 
charge analysis to discriminate between and Ve events, the NOMAD experiment is 
able to limit the systematic uncertainty on the overall near detector flux to 2%, and 



the energy spectrum shape to 2.5% [26|. With a fine-grained near detector of this type 
one might imagine achieving a systematic uncertainty on the far detector flux of 3 to 
4%. 

Background rejection in the far detector can be understood using a second near 
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detector of the same type as the far detector, Assuming that appropriate near detectors 
are used, the remaining systematic uncertainties on the backgrounds in the far detector 
come from the sHghtly different beam spectrum that any near detector sees, and the 
uncertainties on the differences in fiducial volumes and event acceptances. While there 
is probably no hard limit on the systematic uncertainty that could be achieved, it 
is unreasonable to expect that the total systematic uncertainty on the background 
fraction could be reduced below a few per cent of the background itself. In the following 
we will assume that the uncertainty on the background rates is 5%. Note that the 
participants of a recent study of 1 GeV neutrino beams at the JHF @ assumed the 
more conservative value of 10% for the systematic uncertainty on the backgrounds. 

5.6 Summary: Dataset size D & background fraction Jb 

In the previous discussion we concluded that, for a i^^ ^ i^e measurement, the contri- 
bution to the expected background fraction fs from the initial i/g beam contamination 
might be reduced to ~ 0.002. However, for most of the detector types we have consid- 
ered, Jb is dominated by the contributions from (a) neutral pions faking an electron 
signature, and/or (b) i/^ related backgrounds. We would like to know, as a func- 

tion of detector choice, the values of /b and D that should be used in assessing the 

— > physics potential. The dominant contributions to /s, together with fs and 
the associated signal efficiencies, and detector cost estimates, are summarized in Table 
^ along with the implied value of D for 5 years of data taking. The values of D are 
estimated assuming the detectors cost $500M, which determines the detector masses. 
Estimating the unit costs for each detector type is not straightforward. Details of the 
cost estimates are given in Appendix 2. 

Note that no detector achieves the goal fs < 0.004 that we derived in Section 4 
for a multi~GeV beam. Once above the Ur CC threshold the contribution to Jb from 
fsiT) already exceeds 0.004. As an example, consider the liquid Argon detector. The 
parameters to be used for this detector type with a medium energy MINOS-type beam 
are: D = 170 kt-yrs, Jb = 0.008, and crj^/ Jb = 0.05 (as discussed in Section 5.5). If the 
Ur background contribution could be eliminated /b would be reduced to ~ 0.003, but 
this improvement in background rate would be accompanied by a significant decrease 
in D, and would not be expected to significantly improve the sensitivity to 
oscillations. 

6 Physics with multi— GeV long— basehne beams 

Using the values of D and Jb in Table ^, and assuming o-j^/ Jb = 0.05, we can now 
assess the physics potential for the various detector scenarios we have considered. We 
will begin with the multi-GeV long baseline beams, and consider the minimum value 
of sin^ 2^13 that will yield a,v^^V(. signal 3fT above the background, the sensitivity to 
the neutrino mass hierarchy, and the sensitivity to CP violation in the lepton sector. 

6.1 sin^ 2^13 Reach 

To obtain the sin^ 2^13 reach for the detector scenarios listed in Table ^ we return 
to Fig. |5| which shows contours of constant reach in the (/b, -D)-plane for upgraded 
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Table 4: Detector background rates (/s), signal efficiencies, and unit costs. Water cerenkov 
backgrounds and efficiencies are neutrino energy dependent: numbers left of the arrows for 
a 1 GeV beam, numbers right of the arrows for a multi-GeV beam requiring y < 0.5. 





Water Cerenkov 


Liquid Argon 


Steel+readout 


Liquid-^ 




(UNO) 


(ICARUS) 


(MINOS) 


(THESEUS) 


Scintillator 


Signal Efficiency 


0.7^ 0.5 


0.90 


0.33 


o.6^<7 


0.76 


/b(NC) 


0.02 ^ 0.04 


0.001 


0.01 


0.01 


< 0.006 


/^(beam) 


0.002 


0.002 


0.002 


0.002 


0.002 


fB{r) 


^ 0.01 


~ 0.005^ 


0.02'= 


0^ 


~ 0.005^= 


Ib 


0.02 ^ 0.05 


~ 0.008 


0.03 


0.01 


~ 0.01 


Electron cut 


> 0.5 X 


none*^ 


1-6 GeV 


> 0-5Eyis 


E^is > 2 GeV 


Unit cost (M$/kt)'' 


2.4 


23 


10.4 


78 


59 


Mass (kt) / 500 M$ 


745 


37 


85 


6.4 


260 


D (kt-yrs)'' 


2600 ^ 1860 


170 


140 


19 


990 



FYOO dollars. Costs account for salaries, overheads, and contingencies. Details are given 
in Appendix 2. The cost does not include excavating a cavern, which is estimated to be 
0.5M$/kton/p, {p = target density). 
^ For the MINOS low energy beam. 
^ For the MINOS medium energy beam. 

Although a total energy cut might be applied to reduce the intrinsic z/e background. 

For 5 years running. 
^ UeU e~p search, Ahrens et al, Phys. Rev. D volume 36, (702) 1987. 
^ Soudan efficiency for electron neutrinos: NuMI-L-562. 
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Figure 15: Three-sigma error ellipses in the (A^(e"^), A^(e~))-plane, shown for —>■ v^, 
and Vp, — » z/e oscillations in a NuMI-like high energy neutrino beam driven by a 1.6 MW 
proton driver. The calculation assumes a liquid argon detector with the parameters listed in 
Table ^, a baseline of 2900 km, and 3 years of running with neutrinos, 6 years running with 
antineutrinos. Curves are shown for different CP-phases 8 (as labelled), and for both signs 
of 5m32 with 1(5^321 = 0.0035 eV^, and the sub-leading scale ^vn^^ = 10""^ eV^. Note that 
sin^ 26*13 varies along the curves from 0.0001 to 0.01, as indicated [p!3| . 

(1.6 MW) NuMI medium- and high-energy beams. Of the scenarios we have con- 
sidered, the greatest sensitivity is obtained using a liquid argon detector with either 
the medium- or high-energy beams at L = 2900 km, 4000 km, or 7300 km. In these 
cases a Ufj, ^ Ue signal at least 3cr above the background would be expected provided 
sin^ 2^13 > 0.002 to 0.003. If the ly^ backgrounds could be elliminated, reducing /b 
to 0.003, the limiting sensitivity improves to sin^ 2^13 > 0.001. If the initial u^, con- 
tamination in the beam is 0.5% (rather than the assumed 0.2%), so that Jb = 0.01, 
the sin^ 2^13 reach is still ~ 0.002 to 0.003. Hence, the estimated reach is not very 
sensitive to the uncertainties on our background estimations. However, the z/g 
sensitivity would be degraded if the Vr backgrounds were significantly larger, which 
disfavors higher beam energies. 

We conclude that, with a 30-40 kt liquid argon detector and a medium energy 
superbeam, we could improve the sensitivity to — > oscillations, and obtain about 
an order of magnitude improvement in the sin^ 2^13 reach beyond that expected for the 
currently approved next generation experiments. The other detector choices in Table ^ 
do not seem to be competitive with liquid argon, with the possible exception of the 
water cerenkov detector which obtains a reach of about 0.003 with the high energy 
beam at the longest baseline (L = 7300 km). 
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Figure 16: Three-sigma error ellipses in the (iV(/i+), iV(yU— ))-plane, shown for a 20 GeV 
neutrino factory delivering 3.6 x 10^^ useful muon decays and 1.8 x 10^^ antimuon decays, 
with a 50 kt detector at L = 7300 km, 5m2i = 10~^ eV^, and 6 = 0. Curves are shown for 
both signs of ^m^g; sin^ 2^^i3 varies along the curves from 0.0001 to 0.01, as indicated [|T^. 



6.2 Matter Effects and CP Violation 

Having either established or excluded a — > z^e signal, a search for i^^ — > over a 
long baseline can determine the sign of hence the neutrino mass hierarchy. 

Suppose that and A^-_ signal events are measured in neutrino and antineutrino 
running respectively. In the absence of intrinsic CP-violation, and in the absence of 
matter effects, after correcting for different beam fluxes, experimental livetimes, and 
the neutrino/antineutrino cross section ratio, we would expect iV+ = A^_. However, if 
5m32 > 0, matter effects can reduce and enhance A^_. Alternatively, if (Jmgg < 0, 
matter effects can reduce iV_ and enhance A^+. For detector scenarios and multi~GeV 



beams similar to those considered in this report, it has been shown |13] that: 

(i) At L = 732 km the expected changes of and due to matter effects are mod- 

est, and are comparable to changes that might arise with maximal CP-violation 
in the lepton sector. It is therefore difficult to observe and disentangle matter 
from CP effects unless the baseline is longer. 

(ii) At baselines of ~ 3000 km or greater matter effects are much larger than CP 
effects, and the determination of the sign of (5m32 is straightforward provided 
both A'^^ and A_ have been measured with comparable sensitivities, and at least 
one of them is non-zero. The sign of 6m^2 can be determined with a significance 
of at least 3a provided sin^ 2^13 > 0.01. 

(iii) CP violation cannot be unambiguously established in any of the long baseline 
scenarios considered. 

To illustrate these points we choose one of the most favorable scenarios studied: a 
1.6 MW NuMI-like high energy beam with L = 2900 km, detector parameters Jb and 
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D corresponding to the liquid argon scenario in Table ^, and oscillation parameters 
|(5m32| = 3.5 X 10^'^ eV^ and Sm2i = 1 x 10~^ eV^. The calculated three-sigma error 
ellipses in the (A^(e^), A^(e~))-plane are shown in Fig. |l5|for both signs of with 
the curves corresponding to various CP-phases 6 (as labelled). The magnitude of the 
oscillation amplitude parameter sin^ 2^13 varies along each curve, as indicated. 
The two groups of curves, which correspond to the two signs of (5m32, are separated by 
more than 3a provided sin^ 2^13 ^ 0.01. Hence the mass heirarchy can be determined 
provided the — > I'e oscillation amplitude is not less than an order of magnitude 
below the currently excluded region. Unfortunately, within each group of curves, the 
CP-conserving predictions are separated from the maximal CP-violating predictions 
by at most 3cr. Hence, it will be difficult to conclusively establish CP violation in this 
scenario. 

Note for comparison that a very long baseline experiment at a neutrino factory 
would be able to observe — > oscillations and determine the sign of 5m32 for 
values of sin^ 2^13 as small as 0(0.0001) ! This is illustrated in Fig. IC. 



7 Physics with 1 GeV medium basehne beams 

We next turn our attention to neutrino beams with energy E'jy ~ 1 GeV. The atmo- 
spheric neutrino deficit scale 5m'^2 then sets a baseline requirement L ^ 300 km. A 
recent study Q has generated a letter of interest for a 1 GeV neutrino beam at the 
Japan Hadron Facility (0.77 MW 50 GeV proton driver), with a baseline of 295 km, 
using the SuperK detector |31|. The JHF study group has considered a variety of 
beamline designs, including both narrow band and wide band beams, quadrupole and 



horn based focusing. The energy distributions for these beams are shown in Fig. 17 
(left panel). With a water cerenkov detector the JHF group finds that for a - 



e 



search it is important to use a narrow band beam to avoid the high energy neutrino 
tail which provides a significant source of high energy NC events with detected neu- 
tral pions that fake lower energy electrons. With a narrow band beam, after detailed 
studies of signal efficiency and background rejection they obtain fs = 0.03, dominated 
by the surviving NC backgrounds. This background level was obtained at the cost of 
reducing D hy a factor of 0.68. The initial i^e component in the beam contributes only 
0.004 to Jb- An uncertainty on the background rates of 10% was assumed. 

In the following we discuss the sensitivity of a 1 GeV medium baseline neutrino 
beam to — > i^e oscillations, and the prospects for observing CP-violation. Note that 
baselines of a few hundred kilometers are too short for matter effects to be appreciable, 
and hence the pattern of neutrino masses cannot be determined with these medium 
baseline beams. 

7.1 sin^ 2^13 reach 

The JHF study group concluded that, if no fg appearance signal was observed after 
5 years of running, taking the central SuperK value for (J/n-gj, they would be able to 
exclude Ve oscillation amplitudes greater than 0.01 at 90% C.L., which is about 



an order of magnitude better than the present limit (Fig. 17 right panel). This level of 
sensitivity corresponds to a sin^ 2^13 reach of ~ 0.03. With a larger dataset (20 times 
SuperK) the reach is improved by about a factor of 4. 
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Figure 17: Differential event rates (left) and two-flavor v^. oscillation sensitivity (right) 
for a 1 GeV neutrino beam at the JHF with L = 295 km, and a water cerenkov detector. In 
the left panel the thick solid histogram is for a wide band beam, and the thin solid, dashed, 
dotted, and dot-dash histograms are for four narrow band beam designs. The 90% CL. 
contours in the right panel are for, as indicated, 1 year with a wide band beam, 5 years with 
a narrow band beam and the SuperK detector, and 5 years with a narrow band beam and 
20 X the SuperK detector M. 
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Figure 18: Charged pion production for 16 GeV and 50 GeV protons incident on a 30 cm 
copper target. The Tr"*" and vr" kinetic-energy (left) and transverse momentum (right) spec- 
tra are shown for forward going particles within a cone of angle 1 radian. To facilitate a 
comparison at equal beam power the rates have been divided by the incident proton energy. 
For clarity, the 7i~ spectra have been scaled down by an additional factor of 10. 



29 



It is interesting to see where the JHF scenario lies in the (/s, -D)-plane shown in 
Fig. ^. If we take fs = 0.03 and D = 77 kt-years (5 years of data taking in the superK 
detector with a signal efficiency of 0.68), we see that the JHF — > SuperK scenario 
lies close to the sin^ 2^13 = 0.03 contour. Thus the calculations of Ref. [l^ are in 
agreement with the JHF study group results. Note that the JHF scenario already 
lies in the background systematics dominated (vertical contour) region of the plane. 
Upgrading the detector mass by a large factor only results in a modest improvement in 
the sin^ 2^13 reach. With D = 2600 kt-years and /b = 0.02, the reach has improved to 
sin^ 26*13 ~ 0.01 at the 0.77 MW JHF. It is unclear whether an upgraded 4 MW JHF 
would further improve the reach, which is very sensitive to (yj^/fs in the systematics 
dominated region of the (/s, D) -plane. A liquid argon detector, with fs = 0.003, 
Cfg/ /b = 0.05, and D = 170 kt-years, would obtain a reach of ~ 0.01 at the 0.77 MW 
JHF, and ~ 0.003 at an upgraded 4 MW JHF. 

To a good approximation we would expect the JHF study results to apply also to 
a 1 GeV neutrino beam generated at Fermilab using a 16 GeV ~ 1 MW proton driver. 
Charged pion production spectra for 16 GeV and 50 GeV protons are compared in 
Fig. 18, with the spectra normalized by dividing by the proton beam energies. Hence 
the pion event rates are shown at equal beam powers. The shapes of the 16 GeV and 
50 GeV kinetic energy- and transverse-momentum-distributions are similar. At equal 
beam power, the sub-GeV pion rates obtained with 16 GeV protons are approaching 
a factor of two higher than the 50 GeV rates. Above 1 GeV the rates obtained with 
16 GeV protons are similar to those obtained with 50 GeV protons if the beam powers 
are similar. Hence we would expect the 1 GeV neutrino beam fluxes at the JHF to be 
similar to the fluxes at an ~ 1 MW 16 GeV machine. 



7.2 Searching for CP-violation 

If sin^ 2^13 lies within an order of magnitude of the present experimental limits, we 
would expect a f ^ — > signal to be established in either the next generation of 
accelerator based neutrino experiments, or at a future superbeam experiment. In this 
case, if the solar neutrino deficit is correctly described by the LMA MSW solution 
there is the tantalizing possibility of observing CP-violation in the lepton sector, and 
measuring the CP-violating amplitude. In a medium baseline experiment the CP- 
violating signature (an asymmetry between the — > Ue and z^^ oscillation 
probabilities) is not complicated by matter effects, which are very small. 

Consider the sensitivity to CP-violation at the JHF with L = 295 km. Both the 
background levels and the associated systematic uncertainty are expected to be worse 
for antineutrino beams than for neutrino beams. In the following, for simplicity we 
will consider backgrounds and systematics to be the same for v and v beams, and take 
/s = 0.02 and CTfg/fB = 0.1. Figure 19 shows the expected sensitivity to maximal CP- 
violation {5 = 90°) after 3 years of neutrino running to measure the number of — > 
events, followed by 6 years of antineutrino running to measure the number of — > 
events. In the absence of CP-violation we would expect the two signal samples to have 
the same number of events on average since the factor of two difference in neutrino 
and antineutrino cross-sections is compensated by the difference in the running times. 
Hence the broken curves at 45° in the figure correspond to the CP-conserving case. The 
figure shows 3(T error ellipses for a water cerenkov detector with a fiducial mass of 220 kt 
at the 0.77 MW JHF (left panel), and a liquid argon TPC with a fiducial mass of 30 kt 
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Figure 19: Three-sigma error ellipses in the (A^+, A^_)-plane, where is the number of 
i^fj. signal events and is the number of v^i_ —>■ Ue signal events, shown for 1 GeV 

neutrino beams with L = 295 km at the 0.77 MW JHF using a 220 kt water cerenkov 
detector (left panel) with /b = 0.02, and at the 4 MW upgraded JHF with a 30 kt liquid 
argon detector (right panel) with fs = 0.004. The 3 families of ellipses correspond to 
sin^ 26^13 = 0.02, 0.05, and 0.1, as labelled. The solid (dashed) [dotted] curves correspond to 
5 = 0° (90°) [-90°] with Sml^ varying from 2 x 10"^ eV^ to 2 x 10"^ eV^. The error ellipses 
are shown on each curve for = 5 x 10~^, 10""^ and 2 x 10"'* eV^. The curves assume 3 



years of neutrino running followed by 6 years of antineutrino running [|l^]. 



at a 4 MW upgraded JHF (right panel) . The error ellipses are shown for three different 
sub-leading scales 5m2i = 2 x 10^^, 1 x 10^^, and 5 x 10~^ eV'^. In each panel the three 
families of ellipses correspond to three values of sin'^ 2^13 . Note that parameter values 
with ellipses that lie entirely above the 6 = line would result in a 3a observation of 
maximal CP-violation. We see that for the water cerenkov scenario, if sin^ 2^13 = 0.1, 
marginally below the currently excluded region, then maximal CP-violation would be 
observed provided (Jm^i is not significantly below Ix 10~^ eV^, which is at the upper end 
of the preferred solar neutrino deficit LMA region. The sensitivity is only marginally 
better in the liquid argon scenario. With decreasing sin^ 2^13 the sensitivity slowly 
decreases, with 1 x 10^^ eV^ being the limiting (Jm^i for sin'^ 2^13 ~ 0.02. Hence, for a 
small region of presently favored MSW LMA parameter space, maximal CP Violation 
could be seen at 3(T at a 1 GeV medium baseline superbeam. This small exciting piece 
of parameter space can be described approximately by: 

sin^ 26*13 > 0.02 , sin 5 x 5ml^ > 7 x 10"^ eV^ (5) 
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8 Conclusions 



Neutrino superbeams that exploit MW-scale proton drivers, together with detectors 
that are an order of magnitude larger than those presently foreseen, offer the prospect 
of improving the sensitivity to — > Vf, oscillations by an order of magnitude beyond 
the next generation of experiments. Superbeams would therefore provide a useful tool 
en route to a neutrino factory. Our main conclusions are: 

(i) We believe that the initial Ve contamination in the beam might be reduced to 

~ 0.2%, although we note that the contributions from decay will make this 
goal difficult to acheive for multi-GeV beams. 

(ii) The dominant — > t'g backgrounds will arise from (a) vr^ production in NC 
events, where the vr'^ subsequently fakes an electron signature, and (b) v-r CC 
interactions (if the beam energy is above ~ 5 GeV). 

iii) Of the detector technologies we have considered, only the liquid argon detector 
offers the possibility of reducing the background fraction fs significantly below 
0.01. A multi-GeV long baseline superbeam experiment with a liquid argon 
(water cerenkov) detector would be able to observe a — > signal with a 
significance of at least 3cr above the background provided sin^ 2^13 > 0.002 — 0.003 
(0.003). If the baseline L ^ 3000 km, the same experiment would also be able to 
determine the sign of 5'm\2 provided sin^ 2^13 ^ 0.01. However, it seems unlikely 
that an unambiguous signal for CP-violation could be established with a multi- 
GeV superbeam. 

(iv) A 1 GeV medium baseline superbeam experiment with a liquid argon detector 
would be able to observe av^^v^, signal with a significance of at least 3(7 above 
the background provided sin^ 2^13 ^ 0.003. The experiment would not be able 
to determine the sign of 5m\2-, but if the LMA MSW solution correctly describes 
the solar neutrino deficit, there is a small region of allowed parameter space for 
which CP-violation in the lepton sector might be established. 

We compare the superbeam — > reach with the corresponding neutrino factory 
~^ reach in Fig. which shows the 3(7 sensitivity contours in the {dm^i^siv? 2^13)- 
plane. The superbeam sin^ 2^13 reach of a few x 10~^ is almost independent of the sub- 
leading scale (5m2x. However, since the neutrino factory probes oscillation amplitudes 
O(10~^) the sub-leading effects cannot be ignored, and a signal would be observed at 
a neutrino factory over a significant range of dm^i even if sin^ 2^13 = 0. 

Finally, both the possibility of exploiting sub-GeV superbeams (not considered 
in our present study), and the optimum detector designs for GeV and multi-GeV 
experiments, deserve further consideration. 
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Figure 20: Summary of the 3a level sensitivities for the observation of z/^ — > z/g at various 
MW-scale superbeams (as indicated) with liquid argon "A" and water cerenkov "W" detector 
parameters, and the observation of z/g — z/^ in a 50 kt detector at 20, 30, 40, and 50 GeV 
neutrino factories delivering 2 x 10^° muon decays in the beam forming straight section. The 
limiting 3a contours are shown in the (5m|^, sin^ 2^13-plane. All curves correspond to 3 years 
of running. The grey shaded area is already excluded by current experiments. 
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9 Appendix 1: Neutrino Masses and Mixing 



In this appendix we briefly review the theoretical framework used to describe neutrino 
oscillations. 

9.1 Neutrino mass 

In the standard SU(3) x SU(2)lX U(l)y model (SM) neutrinos occur in SU(2)l dou- 
blets with y = -1: 

^Li = (^^^ ^ , £ = e, n, T (6) 

There are no electroweak-singlct neutrinos (often called right-handed neutrinos) XR,ji 
j = l,...,ns. Equivalently, these could be written as X'^lj- There are three types 
of possible Lorentz-invariant bilinear operator products that can be formed from two 
Weyl fermions ipL and XR- 

• Dirac: mo'tpLXR + h.c. This connects opposite-chirality fields and conserves 
fermion number. 

• Left-handed Majorana: mLil^J^Cil^L + h.c. where C = ij2"fo is the charge conju- 
gation matrix. 

• Right-handed Majorana: Mrx^Cxr + h.c. 

The Majorana mass terms connect fermion fields of the same chirality and violate 
fermion number (by two units). Using the anticommutativity of fermion fields and the 
property = —C, it follows that a Majorana mass matrix appearing as 

i;fC{Mmaj)ij^j (7) 

is symmetric in flavor indices: 

MLj = Mmaj (8) 

Thus, in the SM, there is no Dirac neutrino mass term because (i) it is forbidden 
as a bare mass term by the gauge invariance, (ii) it cannot occur, as do the quark 
and charged lepton mass terms, via spontaneous symmetry breaking (SSB) of the 
electroweak (EW) symmetry starting from a Yukawa term because there are no EW- 
singlet neutrinos XR,j- There is also no left-handed Majorana mass term because (i) 
it is forbidden as a bare mass term and (ii) it would require a I = 1, Y = 2 Higgs 
field, but the SM has no such Higgs field. Finally, there is no right-handed Majorana 
mass term because there is no XR,j- The same holds for the minimal supersymmetric 
standard model (MSSM) and the minimal SU(5) grand unified theory (GUT), both for 
the original and supersymmetric versions. 

However, it is easy to add electroweak-singlet neutrinos XR to the SM, MSSM, or 
SU(5) GUT; these are gauge-singlets under the SM gauge group and SU(5), respec- 
tively. Denote these theories as the extended SM, etc. This gives rise to both Dirac 
and Majorana mass terms, the former via Yukawa terms and the latter as bare mass 
terms. 

In the extended SM, MSSM, or SU(5) GUT, one could consider the addition of the 
Xr fields as ad hoc. However, a more complete grand unification is achieved with the 
(SUSY) SO(IO) GUT, since all of the fermions of a given generation fit into a single 
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representation of SO(IO), namely, the 16-dimensional spinor representation ipi. In this 
theory the states XR are not ad hoc additions, but are guaranteed to exist. In terms 
of SU(5) representations (recall, SO(IO) D SU(5) x U(l)) 

16l = Wl + 5l + 1l (9) 

so for each generation, in addition to the usual 15 Weyl fermions comprising the 10^ 
and 5r, (equivalently 5l) of SU(5), there is also an SU(5)-singlet, Xl (equivalently, 
Xr)- So in SO(IO) GUT, electroweak-singlet neutrinos are guaranteed to occur, with 
number equal to the number of SM generations, inferred to be = 3. Furthermore, the 
generic scale for the coefficients in Mr is expected to be the GUT scale, Mqut ~ 10^^ 
GeV. 

There is an important mechanism, which originally arose in the context of GUT's, 
but is more general, that naturally predicts light neutrinos. This is the seesaw mecha- 
nism [p^. The basic point is that because the Majorana mass term is an electroweak 
singlet, the associated Majorana mass matrix Mji should not be related to the elec- 
troweak mass scale v, and from a top-down point of view, it should be much larger than 
this scale. Denote this generically as m/j. This has the very important consequence that 
when we diagonalize the joint Dirac-Majorana mass matrix, the eigenvalues (masses) 
will be comprised of two different sets: Ug heavy masses, of order rriR, and 3 light 
masses. The largeness of rriR then naturally explains the smallness of the masses (or, 
most conservatively, upper bounds on masses) of the known neutrinos. This appealing 
mechanism also applies in the physical case of three generations and for > 2. 

At a phenomenological level, without further theoretical assumptions, there is a 
large range of values for the light rrii,, since (1) the actual scale of niR is theory- 
dependent, and (2) it is, a priori, not clear what to take for m£) since the known (Dirac) 
masses range over 5 orders of magnitude, from me,mu ~ MeV to mt = 174 GeV, and 
this uncertainty gets squared. However, in the SO(IO) GUT scheme, where one can 
plausibly use m/j ~ for the third-generation neutrino, and mji ~ Mqut ~ 10^^ 
GeV for the scale of masses in the right-handed Majorana mass matrix, one has 

2 

m(i/3) ~ ^ ~ 10"^ eV (10) 
rriR 

which is close to the value rn-(t'3) = 0.06 eV obtained from 6rn^^^ if one assumes a 
hierarchical neutrino mass spectrum with m(t'3) >> m[u2). Thus, the seesaw mech- 
anism not only provides an appealing qualitative explanation of why neutrino masses 
are much smaller than the masses of the other known fermions, but also, with plau- 
sible assumptions, predicts a value for 771(^3) comparable to suggestions from current 
atmospheric neutrino data. 

9.2 Neutrino mixing and oscillations 

The unitary transformation relating the mass eigenstates to the weak eigenstates is as 
follows, 

3 

^f-a = X! ' ^1 = e, ^2 = ^3 = T (11) 

i=l 
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I.e., 



( Ve 








1 






\ -3 ) 



U,,l U^2 Uf,3 \\ ^2 \ (12) 
One possible representation of this 3x3 unitary matrix is 

/ C12C13 S12C13 Si3e"*'^\ 

U = -S12C23 - Ci2S23'Sl3e*'' C12C23 - 512523^136"' 523613 (13) 
\ S12523 - 6126235136"^ "612523 " 5126235136"^ 623613 / 

where Cjj = cos 6ij, Sij = sinOij. Thus, in this framework, the neutrino mixing depends 
on the four angles ^12, ^13, ^23, and 5, and on two independent differences of squared 
masses, Sni^i^ , which is 5m'^2 — f^i^s)"^ ~ m{v2)^ in the favored fit, and Sm'^gi, which 
may be taken to be ^m^j^ = m(i/2)^ — m(z/i)^. Note that these quantities involve both 
magnitude and sign; although in a two-species neutrino oscillation in vacuum the sign 
does not enter, in the three species oscillations relevant here, and including both matter 
effects and CP violation, the signs of the 6m? quantities do enter and can, in principle, 
be measured. 

In the 1980's, most theorists thought that lepton mixing would be hierarchical, i.e. 
the lepton mixing matrix U would differ from the identity by small entries, and these 
would be smaller as one moved further from the diagonal, as is established to be the case 
with quark mixing. This was, indeed, a large part of the appeal of the MSW mechanism: 
it could produce large mixing with small vacuum mixing angles. However, the results 
from the SuperK measurements of atmospheric neutrinos have forced a revision in 
this conventional picture, providing strong evidence for essentially maximal mixing, 
sin^ 2^23 = 1- A challenge to model-builders has thus been to get maximal sin^2^23- 
More recently, the SuperK solar neutrino data favors large sin^ 2^12. Bimaximal mixing 



schemes take ^23 = ^12 = vr/4 and ^13 << 1 [33|. There are no compelling theoretical 
suggestions concerning the magnitude of ^13, and one of the important physics goals 
for neutrino oscillation experiments with conventional beams is to try to measure this 
angle. 

For our later discussion it will be useful to record the formulas for the various 
relevant neutrino oscillation transitions. In the absence of any matter effect, the prob- 
ability that a (relativistic) weak neutrino eigenstate Ua becomes Vb after propagating a 
distance L is 

i>j=l 

+ 4 J2 ^"^(^"M.)sm(^^)cos(^^) (14) 

i>j=l 



where 
and 



Kah,ij = UaiUbiU*jUbj (15) 



Sm^j = m{ui) — m{uj) (16) 

Recall that in vacuum, CPT invariance implies P{v'b i^a) = P{T^a ^b) and hence, 
for b = a, P{i>a ^a) = P{va ^a)- For the CP-transformed reaction Va Vb and 
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the T-reversed reaction i/h — > i^a, the transition probabihties are given by the right- 
hand side of ( p!^ ) with the sign of the imaginary term reversed. In the following we 
will assume CPT invariance, so that CP violation is equivalent to T violation. 
The solar and atmospheric neutrino data indicate that 

(^m|i = dml^i <C Smj^ ^ Sml^ = Sml^^ (17) 

In this case, CP (T) violation effects are very small, so that in vacuum 

P{i^a ^ h) - P{l^a ^ l^h) (18) 

P{yb ^ ya) ^ P{ya ^ ^h) (19) 

In the absence of T violation, the second equality ([l^ would still hold in matter, but 



even in the absence of CP violation, the first equality (18) would not hold. With the 



hierarchy (ll7|), the expressions for the specific oscillation transitions are 
P{u,^v^) = 4|C/33p|C/23psin2(^!^) 

= Sin2 2023cos4^i3sm^(^%^) (20) 



P{l^^.^^e) = 4|i7i3p|C/23|'sin2(^ 

= sin2 2%sin2 023sin2(^l||^) (21) 



With units inserted, one has the general relation 

This makes it clear what the approximate sensitivity of an experiment with a given 
pathlength is to a neutrino oscillation channel involving a given Jm^, for a beam with 
an energy E. 

There can be significant corrections to the one-5m^ oscillation formulas if Sml^i is 
at the upper end of the LMA range, STn^^i ~ 10""^ eV^, if sin'^ 2^13 is sufficiently small. 
In this case, keeping dominant terms and neglecting possible small CP violating terms, 
eq. (21) becomes 



P{u,^u,) = sin2 20i3sin2^23sin2(^!^ 

+ sin2 2^12 cos^^cos^ 023 sin2(^^^^) (23) 

Let us denote the two terms as T\ and T2. As an illustrative example, let us con- 
sider a pathlength L sufficiently short that matter effects are not too important. As- 
sume sin^ 2^13 = 0.01 and the upper end of the LMA solution, with sin^ 2^12 = 0.8 
and (5m^Q; = 10^^ eV^. Then for L = 730 km, T2 = O.lTi. For these values, 
s\ii^{5mli^L/{4:E)) = 0.78 while s\t? {Sml^iL / {AE)) = 0.95 x lO^^, so that the path- 
length is causing a strong suppression of the subdominant oscillation due to (Jm^^j. 
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For sufficiently large L and small sin^20i3, the Sra^^i oscillation can be a significant 
contribution to — > Ve- However, we note that making L greater would mean that 
one would also have to make E greater to keep an acceptable event rate with a given 
detector, and this would tend to increase backgrounds to the Va t^e signal. We also 



note that if KamLAND |34| achieves its projected sensitivity, it will have tested the 
LMA solution by ~ 2005. 

In neutrino oscillation searches using reactor antineutrinos, i.e. tests of i>e — > i^e, 
the two-species mixing hypothesis used to fit the data is 

P{U, -^Ue) = l- Sin2 2ereactor sin^ ( '^"^'g"- ^ ) (24) 

where Sm^^g^^-f.^^ is the squared mass difference relevant for z^g — > t^x- In particular, 
in the upper range of values of 5m\^^, since the transitions i>e and Ve 

contribute to Ve disappearance, one has 



P{v, -^v,) = l- sin' 2013 sin^ (^%^) (25) 

i.e., Oreactor = ^13) and the Chooz reactor experiment yields the bound [^] 

sin^ 26*13 < 0.10 (26) 



which is also consistent with conclusions from the SuperK data analysis |31]. 

Further, the quantity "sin^ 20atm." often used to fit the data on atmospheric neutri- 
nos with a simplified two-species mixing hypothesis, is, in the three-generation case, 

Sin2 2eatm = Sin2 2023 COS^ 013 (27) 

The SuperK data implies that (up to redefinitions of quadrants, etc.) 

023 - J (28) 

and sin^ 20i3 << 1. Thus, to good accuracy, 9atm = ^23- 

The types of neutrino oscillations that can be searched for with a conventional 
neutrino beam include: 

• f ^ — > f ^ (disappearance) 

• — > fe, z^e — > e~ (appearance) 

• — > z^T-, z^T- ^ r~ ; T~ — > (e~,/i~)... (appearance) 

Searches for the conjugate oscillation channels require beams. Since these have lower 
fluxes than beams (and this difference can be large with sign-selected vr beams that 
are decaying), one can concentrate on oscillation channels with beams. 

For neutrino oscillation experiments with pathlengths of order 10^ km, matter ef- 
fects are signiflcant. These have been studied in a number of papers, e.g., |3(;, 37, 38 



|39| , 40, 41 1 . The constant density assumption provides a first approximation; realistic 
density profiles were included in the calculations of ^ ^] . In the constant density 
approximation, for a simple two-species mixing, one has 

P{i^a Vb) = sin2(20„) sm'{ujL) (29) 



40 



where 



sin^(26'r, 



sin2(26') + 



UJ 



5m? 



cos(26l) 



V2 



+ 



6m'^ 
~4E 



sm{2e) 



(30) 



(31) 



where is the electron number density of the matter. 



10 Appendix 2: Detector unit costs 

Estimating the unit costs for each detector type is not straightforward. We can base 
our cost estimates on detectors that are currently under construction, or have recently 
been proposed. However, these example detectors have been proposed/costed at differ- 
ent times using different accounting systems in different currencies with different levels 
of external scrutiny. To attempt to compare like-with-like we have started from the 
bare materials and services (M&S) costs of the detector itself, which we have corrected 
to include salaries (SWF), engineering and R&D (EDIA) costs. The scaling factors 
were determined for a current US-based neutrino detector (MINOS). An estimate of 
overheads and contingencies (35%) has been included to reflect the "fully-loaded" costs 
associated with a US-based detector. Finally, the resulting unit costs have been cor- 
rected for inflation to correspond to FYOl dollars. Based on the fully loaded FYOl 
unit costs, for each detector type the mass of the detectors that could be built with a 
budget of $500M can be estimated. The costs for a cavern for each detector technol- 
ogy is based on the recent UNO estimates for a hard rock cavern ($200/m^) using the 
computed detector masses and the densities of the various detector media The 
results are summarized in Table ^. 

The bare detector costs are based on the following: 

(i) The water cherenkov detector estimates are based on those documented in the UNO 

cost estimate [^], and correspond to 237M$/450 kt (FYOO dollars), assuming 
10% photomultiplier coverage in the entire detector. We assume a hard rock site 
rather than the proposed WIPP site. The actual UNO proposal is based on 40% 
photomultiplier coverage in the central third of the detector's volume, which is 
optimized for certain proton decay and astrophysical neutrino channels. These 
physics topics are the driving force behind the UNO proposal. 

(ii) The liquid argon unit cost is based on the Icanoe costs of 14.4 MEuro per 1.9 kt 
module, with $0.9425 per /Euro (FY99 dollars). Note the costing presented in 
this document assumes only cryogenic modules |]28| . 

(iii) The steel-scintillator unit cost is based on the MINOS M&S unit cost, which is 
based on the most recent (2/01) far-detector cost data giving a total of 16.3M$ 



(FY98 dollars) ||, g|. The total far detector (5.4 kt) cost including R&D, labor. 



and institutional overhead costs in then-year dollars is 25.4M$. 

(iv) The mineral oil cerenkov unis cost is based on the MiniBooNE M&S unit cost, 
which is based on TDR detector costs (FYOO dollars) ||l^, 43 1. 
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Table 5: Detector cost estimates. 





Water 


Mineral Oil 


Liquid 


Steel/ 




Cerenkov 


Cerenkov 


Argon 


Scintillator 


Bare Unloaded Unit Cost (M$/kt) 


0.36 


1.3 


7.1 


3.0 


Unloaded Unit Cost (M$/kt) 


0.57 


1.75 


11.2 


4.7 


FY for estimates 


2000 


2000 


1999 


1998 


Loaded Unit Cost (M$/kt) 


0.67 


1.92 


13.5 


5.9 


Mass (kt) per $500M 


745 


261 


37 


85 


Medium density (g/cm"^) 


1.0 


0.9 


1.8 


3.5 


Cavern cost (M$) 


106 


41 


2.9 


3.4 



a) M&S + SWF + EDIA 

b) FYOl costs including overhead and 35% contingency 

c) Note that deep caverns are not necessarily needed. 
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